syms x  a2 R h1;
h1=36.5;
R=21.74;

a2 = asin(x/R) + atan(x/(h1 + R - sqrt(R^2 - x^2)));

ezplot(a2, [0, 17]);
